Defect-induced perturbations of atomic monolayers on solid surfaces 
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We study long-range morphological changes in atomic monolayers on solid substrates induced by 
different types of defects; e.g., by monoatomic steps in the surface, or by the tip of an atomic force 
microscope (AFM), placed at some distance above the substrate. Representing the monolayer in 
terms of a suitably extended Frenkel-Kontorova-type model, we calculate the defect-induced density 
profiles for several possible geometries. In case of an AFM tip, we also determine the extra force 
exerted on the tip due to the tip-induced de-homogenization of the monolayer. 
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Equilibrium properties of monolayers adsorbed on 
ideal, defect-free solid surfaces are by now reasonably well 
understood through a series of experimental and theoret- 
ical works jij. However, most of naturally encountered 
surfaces or surfaces involved in technological processes 
can not be considered as ideal and do contain different 
types of defects, such as, e.g. chemisorbed or adsorbed 
species or surface steps. Experimentally, it has been well- 
known that such defects may have a profound effect both, 
on the adsorption kinetics and on the equilibrium mor- 
phology of the resulting layers. In particular, point de- 
fects often constitute nucleation sites for the adsorbates 
and serve as seeds for island formation Q . On the other 
hand, in the presence of a monoatomic surface step the 
adatoms on the lower terrace are generally attracted to- 
wards the step, which causes their redistribution within 
the layer, as observed, for instance, via intensity oscil- 
lations of thermal He scattering at grazing incidence in 
the form of the discrete row growth of Xe on stepped 
substrates || ||]. Theoretically, the impact of the sur- 
face steps on the adatom distribution was studied within 
the framework of two-dimensional (2D) lattice-gas-type 
models j5| g|. These models have been analysed numer- 
ically and have revealed inhomogeneous density profiles 
with an enhanced density close to the lower step edges. 
To the best of our knowledge, however, the analytical 
solution of the problem is still lacking. 

On the other hand, probing of the monolayer proper- 
ties by different experimental techniques, such as, e.g. 
the STM or AFM measurements, may itself incur mor- 
phological changes into the adlayer. The interaction of 
the adatoms with the AFM tip might cause their dis- 
placement from the adsorption sites. Such deforma- 
tions have been predicted for solid surfaces themselves 
and were indeed observed in MD simulations ||. The 
adatoms are, or course, even more vulnerable to the pres- 
ence of the AFM tip since they are not so strongly con- 
nected as the atoms of the solid. Indeed, it has been 



demonstrated that the SFM tip can be used to "drag" 
single atoms or molecules on metal surfaces J9|, [lC| [ll], [l2| . 
Moreover, it has been observed in recent experiments 113 ] 
that the apparent thickness of the prewetting film on a 
silicon wafer measured by the AFM is larger as the one 
found by an X-ray reflectivity experiment. The authors 
thus concluded that the AFM tip distorts the film and 
induces a "bump" in its surface. At sufficiently high tem- 
peratures, even stronger effects like the formation of a 
neck between the adsorbate and the tip have been ob- 
served Q. Usually, such a distortion of adlayers is un- 
accounted for while interprctating the experimental data, 
although its effect might be not negligible - "condensa- 
tion" of the adlayer particles in the vicinity of the tip 
would increase the force exerted by the monolayer on 
the AFM tip. Thus the question arises of how to in- 
terprete the AFM measurements adequately and how to 
extract, in a reliable fashion, the pertinent parameters 
(say, the Hamaker constants) in case when some adsor- 
bate is present on the solid surface. 

In this paper we study perturbations of atomic mono- 
layers on solid substrates induced by immobile defects 
of different types. The monolayer is described using a 
2D version of the Frenkcl-Kontorova (FK) model, i.e., 
we view it as a 2D network of particles connected by 
harmonic springs in a spatially periodic potential. Note 
that the original FK-model (a harmonic chain in a spa- 
tially periodic potential) was introduced more than sixty 
years ago in order to describe the motion of a disloca- 
tion in a crystal |jl5| . In the meantime variants of this 
model were applied to many different problems includ- 
ing charge density waves sliding friction jl^, [l8| , 
ionic conductors | IE , |2(J and chains of coupled Joseph- 
son junctions (l6[ 21 1 . A 2D version of the FK-model has 
been introduced by Uhler and Schilling to study glassy 
properties of an adsorbed atomic layer |22 |. 

We consider first the case of a surface with a 
monoatomic step focusing on two opposite limits: (a) 
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monolayers with strong intralayer coupling and negligi- 
ble interactions with the substrate (smooth structureless 
surface) and (b) monolayers in which coupling to the sub- 
strate dominates the particle-particle interactions. As a 
second example we calculate the perturbation of a mono- 
layer induced by an immobile AFM tip and demonstrate 
how it modifies the force exerted on the tip. 

Note that considering the simplified case with an im- 
mobile AFM tip allows us to determine explicitly inhomo- 
geneous density profiles as well as to elucidate the physics 
behind this effect. In "real stuff" experiments the tip, of 
course, moves and the situation is fairly more complex 
since the density profiles are non-stationary. One ex- 
pects that for moderate tip velocities the "condensed" 
region in the monolayer would travel together with the 
tip exerting some frictional force on it. For larger veloci- 
ties, such a condensed region would not have enough time 
to be formed and the monolayer should remain homoge- 
neous. One expects hence the existence of a threshold 
tip velocity below which the monolayer has time to re- 
organize itself leading to an extra force and above which 
this effect disappears, a type of force-velocity relation 
that is somewhat reminiscent of "solid friction" behav- 
ior. A qualitatively similar behavior has been predicted 
by Raphael and de Gennes [^3) for a system involving 
a charged particle moving at a constant speed a small 
distance above the surface of an infinitely deep liquid. 
The situation with a stationary moving AFM tip will be 



discussed elsewhere (24 



We constrain ourselves here to the limit of the local- 
ized adsorption pj and suppose that the adatoms al- 
ways remain in close contact to the surface, such that 
their defect-induced displacements in the vertical direc- 
tion (perpendicular to the surface) are negligibly small. 
For simplicity, we assume that the adatoms form a reg- 
ular square lattice. Each particle is labeled by two inte- 
gers, (n,m), with n,m = 0,±f,±2, .... For small pertur- 
bations of the monolayer the interaction between a given 
atom (n, m) and its four neighbors (n — 1 , m), (n + 1, m), 
(n, m — I) and (n, m + 1) can be represented by Hookean 
springs that connect each atom to its neighbors. The 
value of the effective spring constant K follows from 
the expansion of the interaction potential between atoms 
near the equilibrium distance and is typically of the order 
of a few tenth of eV/A 2 [[f8) . In the absence of any exter- 
nal perturbation the position of atom (n, m) is given by 
the two-dimensional vector r nm — (x nm ,y nm ) — (bn,bm) 
with b being the equilibrium distance between atoms. In 
the following we calculate the defect-induced displace- 
ments a„ m = (£„ m ,77„ m ) of the adatoms. 

First we consider the equilibrium properties of a mono- 
layer near a steplike defect (Fig. f). The substrate has 
the height z = for x > and z — h for x < 0. We 
focus on the monolayer on the lower terrace (x > 0). A 
given small volume element dV of the substrate is as- 
sumed to exert a force df = — adVr~^ a+1 ^r on a particle 




FIG. 1: Atomic monolayer at a steplike defect of height h. 
The particles are attracted towards the step and the mono- 
layer is perturbed accordingly (see text for details) 



in the monolayer at a distance r apart; here a is a con- 
stant, a + 1 is an arbitrary positive number and f is the 
unit vector in the r-direction. In the absence of a step, 
h = 0, the interaction of any atom with the substrate 
is isotropic with respect to rotations around the Z-axis. 
Hence, by summing over all forces between a given atom 
and the substrate one finds only a force perpendicular to 
the surface but no tangential component. In this case 
the monolayer is unperturbed. On the other hand, a 
step of height h > results in net forces to the left for 
atoms to the right of the step. Consider an atom at 
r = (x, y, z) — (D > 0, y, 0). The net force follows from 
integration over the additional slab of material: 
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with C a = V^r((a + l)/2)/((a-I)r(I + a/2)), 
where T(z) is the Gamma-function. The right hand side 
of Eq. [I] holds for a > 1 and h <C D, The first condi- 
tion is needed to insure that /( ste P) remains finite. The 
second condition is fulfilled for small step heights {e.g. 
monoatomic steps, h ~ 6). 

The positions of the atoms obey the force balance equa- 
tion x n j r \ ra 2x nm T x n — \ m A ^ (x nra ) that 
can be rewritten using continuous variables n and m 
as: d 2 x n /dn 2 = (C a ah/ A a;" -1 ) . We drop here the in- 
dependence since the problem is apparently symmetric in 
the Y-direction; x n (n = 1, 2, ...) denotes the apposition 
of the n-th row of atoms in the monolayer. We denote by 
£n = x n — bn the displacement of the n-th row. Assuming 
in the following a weak perturbation of the monolayer, 
b, we are led to 
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where I has the dimension of a length. The exact position 
of the first row (and therefore the value of the constant 
B) depends on its specific interaction with the step. One 
may, for simplicity, assume £i = and thus B = —I. 
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Note that due to the " coupling" of the different rows the 
displacement increases with n approaching the limit- 
ing value B = —I. For a = 6 (van der Waals interaction) 
we find from Eq. §: £„ ~ (tt/480) (ah/ Kb 5 ) (n~ 7 - l). 

The density of adatoms follows p n — b~ 1 dn/dx n — 
b~ 2 (l — b~ 1 d£ y n/dn), which leads to 



C a ah 1 
aKb a n a 



(4) 



i.e. a long-range algebraic relaxation to the un- 
perturbed density. For a — 6 we find p n ~ 
b~ 2 (l + (nah) / (96K [bnf\ ) . The density has its max- 
imal value at the step and decreases with increasing n. 

Up to now we assumed only an intralayer interaction 
between atoms in the monolayer. The role of the sup- 
porting substrate was only to restrict the motion of the 
atoms in the XY-plane. Now we study the case of strong 
coupling to the substrate. We assume next that each 
particle in the monolayer is attached to the substrate 
via a spring with the spring constant K at the equilib- 
rium position r nm = (x nm ,y nm ) = (in, bmj (b might be 

considered to be the lattice constant of the substrate). 
We neglect interactions between neighboring beads, i.e., 
we set K = 0. Then the displacement of the particle 
rows is obtained directly from the force balance equation 

K£ n ~ f (step *> (bn\ . We find then 
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The perturbation of a monolayer with strong intralayer 
coupling (and negligible coupling to the substrate) is fun- 
damentally different from the case of strong coupling to 
the substrate. In the first case - due to the connectivity 
of the rows - the displacement of each row adds up and 
the largest displacement, I (in negative .^-direction) , is 
approached for large n- values, i.e., far from the step (cf. 
Eq. ||) . In the latter case the displacement is directly pro- 
portional to the exerted force which decays algebraically 
with increasing distance and — > for n — > oo. This 
is also reflected in the density profile. The density of 
the monolayer with intralayer coupling, Eq. [|, has its 



largest value close to the step and decays towards the 
unperturbed value b~ 2 for large n. On the other hand, 
monolayers coupled to the substrate show a slight de- 
pletion (p n < b~ 2 ) close to the step, cf. Eq. ||. Only 
very close to the step (first row of atoms) the density 
is enhanced accordingly. The general case with non- 
vanishing K and K is highly non-trivial (e.g., in view of 
possible commensurate-incommensurate transitions) and 
is beyond the scope of the present paper. 

We study next the perturbation of the monolayer by an 
AFM-tip located at height H above the "central" atom 
n = m = 0. We assume that the interaction energy 
between the tip and a particle at distance d apart is of 
the form w (d) = —Ad~ a , which yields 
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where d nm = \(x nm , y nm , H)\ and r nm is the 2D unit 
vector (x n ,y n ) / \ (x nm ,y nm )\. 

The calculation of the elastic force in a monolayer with 
intralayer coupling is non-trivial, since the equilibrium 
distance of each spring has a non-vanishing value b > 
0. For b = the elastic force is simply given by the 
Laplacian: flm™"^ = K (d 2 /d 2 n + d 2 /d 2 m) r nm . For 
b > the X and Y direction are coupled in a non-trivial 
way. However, the elastic response to small perturbations 
a nm = (£nm> Vnm) with |a nm | < b decouples in the X and 
Y directions: 
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The particle positions follow from the balance between 
the tip-monolayer interaction, Eq. ^ and the elastic force, 
Eq. 0, which gives 
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where we have introduced the dimensionless parameter 
7 = H/b. Further on, the displacement £ nm has to be 
calculated for each a separately. For a — 6 (van der 
Waals forces), for instance, one finds: 



-A ( 3n 2 + 5 (m 2 + 7 2 ) / n \ 3 \ 

48Kb 7 \ U ( m 2 + 7 2) 2 ( n 2 + m 2 + 7 2) 2 + arCtai1 ^ ( m 2 + 7 2)V2 ) ^= 2 = J 



Note that a m „ is not radial-symmetric around 
(ri = 0, m = 0) even though it is induced by a radial- 
symmetric force, Eq. |[ In fact, for large n £ nm oc n° 



for m = (X-dircction) and t; nm oc n~ 5 for m — n (di- 
agonal direction). Such a non-isotropy appears due to 
the symmetry of the underlying lattice (see Fig. 2). 
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FIG. 2: View from above on a monolayer close to an AFM 
tip. The tip is located at the height H above the central atom 
(n = m = 0); see text for details. 



For small deformations the density profile of the mono- 
layer is given by 
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Note that here, however, despite the asymmetry of a m „, 
the resulting density profile, Eq. ^|, recovers the symme- 
try of the force exerted by the tip. 

We calculate next the force F that the monolayer ex- 
erts on the AFM tip. Due to the symmetry this force 
is pointed into the negative Z-direction. A particle 
at (x,y) contributes to this force by J^ 4 ™' ,p ' (x,y) = 
—A/Hd a+2 . The total force F from the monolayer fol- 
lows by summing up over all atoms, F = F + AF where 
Fq is the force that an unperturbed monolayer would ex- 
ert on the AFM tip, F = -2tt A/ ab 2 H a - 1 , while AF 
denotes the contribution due to the self-induced pertur- 
bation by the monolayer AF = -2nA 2 /a 2 Kb A H 2a - 1 . 
Note that AF/F = -F /2nKH, i.e., the induced addi- 
tional force is important for soft monolayers (small K) 
and strong tip-sample interactions. For the case a — 6 
the two contributions to the force are given by Fq ~ 
- (vr/3) (A/b 2 H 5 ) and AF ~ - (tt/18) (A 2 /Kb 4 H n ). 

We consider now the effect of small surface corruga- 
tions of the form [/( sur /) (x, y) — eUq cos (kx) cos (ky) on 
the positions of the atoms in the monolayer as well as 
on the force on the AFM tip. Here k denotes the wave 
vector of the periodic substrate and e is a small number, 
e <C 1. From the potential follows the force that acts on 
a particle at (x, y): Fj( Ur (x, y) — — <9[/( SMr /) (x, y) jdx. 
We calculate the additional displacement due to the cor- 



rugations using the ansatz 
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As a result of this additional perturbation the force on 
the AFM tip will be modified, F = F<°) + eFW with 



F(°) = F + AF given above. We give here explicitly the 
asymptotic forms of F^ for the case a = 6: 
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for kH < 1 
' V2kH for kH > 1 



(12) 



It can be seen from Eq. [12| that the contribution from 
surface corrugations is "screened" when the height H of 
the AFM tip exceeds the wavelength of the corruga- 
tions. We dispense with giving a discussion of the case of 
strong coupling to the substrate which can be calculated 
straightforwardly. 
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